Cosmic gamma-ray background from Type la supernovae revisited: Evidence for 

missing gamma-rays at MeV 
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The observed cosmic 7-ray background at ~ MeV has often been attributed to Type la supernovae 
(SNIa). Since SNIa is close to a standard candle, one can calculate the 7-ray intensity of SNIa 
integrated over redshifts fairly accurately, once the evolution of the SNIa rate is known. The latest 
SNIa rate measured at z < 1.6 (Dahlen et al. 2004) indicates that the previous calculations of the 
7-ray background consistently overestimated the SNIa rate. With the new rate, we find that the 
SNIa contribution is an order of magnitude smaller than observed, and thus new population(s) of 
sources should be invoked. 

PACS numbers: 95.35.+d, 95.85.Nv, 95.85.Pw 



The cosmic X-ray and 7-ray background encodes the 
most energetic phenomena in the universe. It has widely 
been accepted that different population of sources con- 
tribute to the different energy bands (see for a re- 
view). On the low energy side, the intensity spectrum in 
the X-ray or soft 7-ray region, < 0.5 MeV, has been suc- 
cessfully explained by the integrated counts of obscured 
Active Galactic Nuclei (AGNs) HH0. At ~ °' 5 MeV 
the spectrum of these AGNs sharply cuts off 4]. On 
the high energy side, 30 MeV < E < 10 GeV, beamed 
AGNs (blazars) are able to account for the observed back- 
ground almost entirely 0, 0, 0- The blazar spectrum 
appears to break below ~ 10 MeV; thus, blazars are un- 
able to account for the low energy spectrum. A natu- 
ral question to ask is, then, "what are the most domi- 
nant sources contributing to the medium energy band, 
0.5 MeV < E < 10 MeV?". 

Ever since the first proposal made by Clayton and Silk 
in 1969 8J, it has often been argued that the 7-ray back- 
ground at ~ MeV region can be accounted for by Type 
la supernovae (SNIa) H d EJ El The contribution 
from core-collapse supernovae at ~ MeV must be much 
smaller than that of SNIa because the 7-ray photons can- 
not easily escape the hydrogen envelope of the progenitor 
(a massive star) E3- If this is true, the spectrum of the 
7-ray background can be used as a powerful probe of the 
cosmic star formation history. The previous calculations 
of the 7-ray background from SNIa were, however, sub- 
ject to uncertainty in the supernova rate (SNR) of SNIa. 

In this paper, we present a more robust prediction for 
the cosmic 7-ray background from SNIa using the obser- 
vationally determined SNR, thereby avoiding any uncer- 
tainty associated with the progenitor model. We show 
that the previous calculations consistently overestimated 
the SNR, and the SNIa contributes no more than 10% 
of the observed level of the cosmic 7-ray background. 
Our results strongly argue for the existence of new pop- 
ulation^) of sources contributing to the background at 
- MeV. 

The spectrum we derive in this paper was already pre- 
sented in the previous paper [l3l | where it is shown that 



the soft 7-ray background at < 0.511 MeV may have a 
substantial contribution from the redshifted 0.511 MeV 
lines of dark matter annihilation in galaxies distributed 
over cosmological distances. These signals cannot, how- 
ever, explain missing 7-rays at > 0.511 MeV. 

We calculate the background intensity, I v , as E3 
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where v is an observed frequency, H(z) is the expansion 
rate at redshift z, and P v {y,z) is the volume emissivity 
in units of energy per unit time, per unit frequency and 
per unit proper volume: 



P v [y,z) = {l + zf SNRi a (z)^. 
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Here, SNRi a is the SNR of SNIa, which is the number 
of SNIa per unit time and per unit comoving volume. 
(Hence (1 + z) 3 in the front.) The time-averaged 7-ray 
energy spectrum of each supernova, E v , in units of en- 
ergy per unit frequency, depends only on v as SNIa is a 
standard candle to a good approximation. One obtains 
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where £ 2 {z) = fl m h 2 (l + z) 3 + Q\h 2 . It follows from 
equation © that contribution from high redshift SNIa 
is negligible unless SNRi a (z) grows faster than (1 + z) 5 / 2 
toward higher redshifts. On the other hand, recent obser- 
vations suggest that SNRi a (z) peaks at z ~ 1 and drops 
at higher redshifts; thus, the most dominant contribution 
must come from lower redshifts. Dahlen et al.|l5j have 
determined the SNRi a at z < 1.6 based on 25 SNIa ob- 
served as part of the Great Observatories Origins Deep 
Survey (GOODS). Figure ^ shows the data (symbols) as 
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FIG. 1: Type la supernova rates (SNRi a ). Observed data 
from GOODS survey are shown with la statistical error bars 
(filled circles). The solid line is the best-fit model deduced 
by Strolger et al. ll6H. The dashed line is for the rates used 
in Watanabe et al.[l2l|. The model adopted by The et al.[ld| 
gives the SNIa rate of ~ 1.6 x 10 -3 Mpc -3 yr -1 , which is 
an order of magnitude larger than the best-fit model and 
is not shown. The dotted line is constructed by connecting 
the 2-cr upper limits (including systematic errors) of observed 
rates[l5l|. 



well as the best-fitting model of Strolger et al.[T^ (the 
solid line). The SNRi a (z) reaches 1.6 x 1CT 4 Mpc" 3 yr" 1 
at 0.6 < z < 1.0 and then drops at z > 1. It thus fol- 
lows from equation J2J) that the expected intensity of the 
7-ray background is vl v ~ 0.5 keV cm~ 2 s _1 str" 1 , as 
£(z)(l + z) ~ 2 at z = 0.8. 

As a comparison, we also plot in Figure^the SNRi a (z) 
used by the previous work|12j. The latest determina- 
tion by lies below it essentially because their es- 
timate was based upon indirect determinations. They 
used an empirical analytic function for the SNRi a (z) 
normalized to the local determination of the star for- 
mation rate at z = 0, assuming a Type II supernova 
rate per unit stellar mass of 0.007 M^ 1 and Type la to 
Type II ratio of 1/3. Another uncertainty comes from 
an adopted progenitor model of SNIa. Observations sug- 
gest that SNRi a (z) follows the star formation rate with 
a delay time. This is where significant uncertainty in 
the progenitor model comes in. They showed that the 
uncertainty in the delay time easily affects the ampli- 
tude of the predicted 7-ray background by a factor of 2 
to 3 (Figure 5 of ^3)- (The longer delay time reduces 
the amplitude.) The SNRi a used by [lfj is even larger 
(SNRi a ~ 1.6 x 10~ 3 Mpc~ 3 yr -1 ) and is not shown in 
Figure n They used an analytic model of the SNRi a (z) 
for which it was assumed that all the baryons in the uni- 
verse were converted into stars. This led to a gross over- 
estimation of the star formation rate by an order of mag- 
nitude. Now that SNRi a (z) has been determined up to 
z ~ 1.6, we can circumvent all of these theoretical un- 
certainties in the star formation rate and the progenitor 



model, by taking observationally determined SNRi a (z). 

Currently, the most favored scenario for SNIa is the 
explosion of a single degenerate White Dwarf close to 
the Chandrasekhar mass, in which the ignition occurs 
close to the center. Initially, the nuclear burning front 
propagates as a deflagration front and, after burning of 
about 0.2 to 0.3 M Q , turns into a detonation 01 ■ This 
so called "delayed-detonation" model allows to explain 
both the optical and infrared light-curves and spectra. 
This class of models produce the wide range in masses of 
radioactive 56 Ni between 0.08 to « 0.8Af Q needed to ex- 
plain both "normal-bright" and subluminous SNIa, and 
the brightness decline-relation 0] , a cornerstone of mod- 
ern cosmology. For recent reviews, see 0, |2jj. As a 
typical feature of these models (and consistent with late 
time spectra |2l|), the mean 56 Ni velocities increase with 
the 56 Ni mass and, thus, the relative contribution of SNIa 
with high 56 Ni masses dominate the X-ray spectra, avoid- 
ing potential problems caused by the uncertainties in the 
rate of subluminous SNIa j3^. By using a SNIa at the 
bright-end of the distribution, we may overestimate the 
SNIa contribution to X-rays. 

As a reference model for the 7-ray spectrum and 
the light curve of SNIa, we have chosen the delayed- 
detonation model, 5p0z22.23, of which produces 
0.561M© of 56 Ni and represents a "typical" SNIa, i.e., 
at the bright end, where most SNIa are observed [39| . 
The calculations are based on a Monte Carlo code for 
7-rays [52] but with updated bound-free opacities and 
nuclear branching ratios j2|| [53] • 

Figure |2 shows the time-averaged 7-ray spectrum of 
SNIa. The lines correspond to various 7-ray emission 
lines from radioactive decays of 56 Ni — > 56 Co and 56 Co — > 
56 Fe. The line energies and branching ratios are summa- 
rized in Table 1 and 2 of [H- All lines above 1 MeV 
except one at 1.5618 MeV are due to the decay of 56 Co, 
whereas all lines below 1 MeV except one at 0.84678 MeV 
are due to 56 Ni. 

Figure [3] shows the predicted spectrum of the cosmic 
7-ray background from SNIa as well as the observational 
data from HEAO-1 A4 MED^ and COMPTEL0 ex- 
periments. The AGN contribution 4], which explains the 
HEAO-1 data at lower energy, is also shown. It is quite 
clear that the predicted SNIa signal falls short of the 
COMPTEL data by at least an order of magnitude. We 
also plot the 2-cr upper limit of SNRi a (z) [l5j (includ- 
ing systematic errors), finding that the 2-cr limit is still 
a factor of 6 smaller than observed. Therefore, we con- 
clude that SNIa cannot account for the observed 7-ray 
background at ~ MeV, and other sources should be in- 
voked. Similar results were obtained independently by 
plj , using their "concordance star formation rate" . 

To the best of our knowledge, there are no confirmed 
sources which could produce a substantial amount of the 
cosmic 7-ray background in this energy region. We thus 
argue that there should be new population(s) of sources 
accounting for the "missing 7-rays" at ~ MeV. What 
could these sources be? Perhaps the most straightfor- 
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FIG. 2: top: Time-averaged supernova spectrum, E u , of a 
proto-type SNIa. bottom: The same spectrum expressed in 
differential photon number (=E v /{hv)). 




log |0 MkeV) 

FIG. 3: Predicted 7-ray background from Type la super- 
novae. The solid line is the prediction from the best-fit SNRi a 
by |T3 |. while the dotted line is the one using the 2-a upper 
limit (including systematic errors; see also Figure 0. The 
data points with error bars are from the HEAO-1 A4 MED|2S[ 
and COMPTEL experiments^. The AGN contribution^, 
which explains the HEAO data, is also plotted in the thin 
dashed line. The sum of the SNIa and AGN contributions is 
plotted in the thick lines. The top and the bottom panels use 
different units. 



ward possibility would be a population of blazars emit- 
ting in ~ MeV region. These "MeV blazars" [3(| cannot, 
however, be the primary candidate; otherwise one must 
require all of the regular blazars to be the MeV blazars, 
contrary to observations [3l| . 

If one had to abandon "ordinary" astronomical sources 
such as AGNs and SNIa (which we argue we should) as 
an explanation to the MeV 7-ray background, then more 
exotic sources would be required. Of potential candidates 
would be 7-rays from dark matter annihilation, although 
popular dark matter candidates (e.g., neutralinos) are 
usually very heavy (dark matter mass of m x > 30 GeV) 
and it is unlikely that such heavy dark matter parti- 
cles contribute to the MeV 7-ray background 32] . Much 
lighter dark matter (1 MeV < m x < 100 MeV), on the 
other hand, is more promising, and it has recently been 
shown |33j| that such light dark matter is a promising ex- 
planation to the 511 keV lines detected at the center of 
our Galaxy jsU, |3^] . We have shown in the previous pa- 
per that the redshifted 0.511 MeV lines from other galax- 
ies distributed over cosmological distances contribute to 
the cosmic 7-ray background at < 0.511 MeV substan- 
tially [l^ • Although these lines do not produce any flux 
at > 0.511 MeV, it may not be so surprising that any 
other associated continuum emission, such as the internal 
bremsstrahlung H^, might also contribute substantially 
at ~ MeV Although much is uncertain, the MeV 

7-ray background would certainly be a potential window 
to the new physics. 

Finally, we emphasize the fact that the precise shape of 
the spectrum of the MeV 7-ray background is currently 
not very well constrained (see Figure [SJ. Given the im- 
portance of this region of the spectrum, it seems urgent 
to carry out more precise measurements of the MeV 7-ray 
background. 
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